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A simulation study using the time-dependent Ginzburg-Landau theory is performed for the vortex 
state in two-component superconductors, such as PrOs4Sbi2. We investigate the flux flow and the 
pinning of the vortex sheet structure. We find domain wall that traps half flux-quantum vortices 
and moves with the flux flow. In the pinning case, we observe an emitting process of a conventional 
vortex from the vortex sheet by combining a pair of half flux-quantum vortices. 
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Some exotic unconventional superconductors have 
multiple superconducting phases, indicating a multi- 
component pair potential of the superconductivity. In 
these superconductors, we expect interesting behavior to 
occur through a combination of the multiple components. 

In the newly discovered superconductor PrOs4Sbi2, 
which is a heavy fermion compound with a filled skut- 
terdite structure, thermal conductivity experiments have 
shown a multiple superconducting phase diagram with 
two phases, the A (in the higher H region) and B (lower 
H) phases.^ Multi-component superconductivity is also 
suggested by the spontaneous moment in the supercon- 
ducting phase reported in muon spin relaxation (/iSR) 
experiments.^ As an another example, UPts has a dou- 
ble transition of superconductivity and the superconduct- 
ing phase is divided into three phases, indicating multi- 
component superconductivity. ^ '^'^ 

Considering the two component case for simplicity, the 
pair potential in a 2 x 2 matrix form can be decomposed 
as A(r, k) = 77i(r)(/)i(k) + 7^2(1*) 02 (k) with the order pa- 
rameter 7?m(r), where r is the center of mass coordinate 
of the Cooper pair and m = 1,2. The pairing function, 
depending on the relative momentum k of the pair, is 
given by (j)mO^) = io-y(j)rnO^) ^^r the spin singlet pair- 
ing, and ^m(k) = i^j=x^y^z^rnjO^)o-jd-y for the triplet 
pairing with Pauli matrices , , and . In the double 
transition of superconductivity, after the first component 
?7i appears at the first transition, the second component 
r]2 appears at the second transition. We discuss in this 
paper the exotic vortex structure due to the combination 
of the two components below the second transition. 

Superfluidity of ^He, which has A and B phases, 
is also a typical example of multi-component pair po- 
tentials. Chiral p-wave superconductivity in Sr2Ru04 
should be treated as a two-component superconductor 
with kx + iky-wdive and kx — iky-waye components, or 
kx-wdive and ky-waxe components, when we consider the 
vortex, ^^'^^-^^ the domain wall and boundary states. 

In multi-component superconductors, some supercon- 
ducting states are degenerate in free energy. In a simple 
example of weak coupling singlet pairing with = <pA 
and (j)2 = i^B, two states A± = t^i^a ± ^^20b are degen- 
erate when |A+| = |A_|. Therefore, a domain structure 



may appear in a multi-component superconductor, i.e., 
some regions in a sample are a domain and others 
are a A_ domain. Between the domains, domain walls 
appear as topological defects, which are not easily de- 
stroyed. 

When magnetic fields are applied to this domain struc- 
ture, some of the vortices are trapped at the domain wall. 
The vortices at the domain wall form an exotic structure 
called a "vortex sheet", where a conventional vortex 
splits into two vortices with half flux- quant a ^SiiSiii We 
confirmed the appearance of the vortex sheet by a simu- 
lation in a previous paper. ^ This interesting vortex sheet 
structure can be a clue to the presence of the domain 
wall, which would be clear evidence of unconventional 
multi-component superconductors. 

When a super current flows in the vortex state, vortices 
flow in the direction transverse to the current, which is 
the origin of the flux flow resistivity. In order to stop 
the vortex motion and prevent flux flow resistivity, we 
need to introduce pining centers for vortices in the sam- 
ple. The purpose of this paper is to investigate the dy- 
namics of the vortex sheet structure by simulation of flux 
flow and pinning based on the time-dependent Ginzburg- 
Landau (TDGL) theory, and to examine how vortices and 
the domain wall of the vortex sheet structure move un- 
der the current flow in the presence of pinning centers. 
These dynamics may affect the flux flow resistivity or the 
magnetization process in multi-component superconduc- 
tors. The TDGL theory is used as a phenomenological 
approach in our qualitative study, with the expectation 
that vortices move so as to approach the free energy min- 
imum state. Contrary to a full-gap superconductor, since 
unconventional superconductors have low energy quasi- 
particle states available for the dissipation process within 
the superconducting gap due to the gap node, we expect 
that the TDGL theory can be qualitatively applied with- 
out considering the extreme case of the dirty limit. 

We start from the two-component Ginzburg-Landau 
(GL) theory that was used in our previous study on the 
double transition of PrOs4Sbi2. The GL free energy den- 
sity in the dimensionless form is written as^- 
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^Ci2xViqIv2 + C^2xV2qIvi 

^Ci2yVUlV2 + Cl2yV2QlVl 

+ ^{|^i|' + C'2|r?2|^ + 4C3|r/ip|r?2p 

+C4*7?f7?i+^4/?2*'^?}, (1) 

where q = (Ti/i)V + (27r/^o)A with flux-quantum 
and vector potential A, and we write C22a; = (1 — 
e)/{X^/T^) and ^222; = {I ^ e) / {X Vl^) . Inside 
the pinning center region, where the superconductivity is 
suppressed, we set z^(r) = 0.^^ In other regions, v{y) = 1. 
The GL theory covers both singlet and triplet pairings 
and we can study the vortex states without specifying 
the pairing function ^(k). The coeflicients in Eq. Q 
are related to the information of the pairing symmetry 
and the Fermi surface anisotropy as described in Ref. 9. 
However, since this information has not been established 
yet, we treat the coefficients as arbitrary parameters. 

We consider a phase below the second transition, where 
both 771 and 772 appear. From Eq. Q, we recognize that 
the relative phases of r]i and r]2 should be (a + 7r)/2 or 
(a — 7r)/2 in the free energy minimum state in a zero 
fleld, where a is given by C4 = |C4|e*". These two states 
correspond to A+ and A_ noted in the introduction. 

For simulations of the vortex dynamics we use the 
TDGL equation coupled with Maxwell equations 

1 1 (2) 

dt^^ 12 St?*' dt^^ 12^7?*' ^ ' 

^A=j,-K^V xB, B = VxA. (3) 
ot 

The supercurrent js = {js,xjs,y) oc {df/dA^^df/dAy) is 
given by 

js,x = Re[r]l{qxr]i) + C22xr]2{QxV2) 

-^Ci2xVi{QxV2) + C^2xV2{QxVi)]^ (4) 

js,y = R^ViiQyVi) + C22yV2{QyV2) 

-^Ci2yVUQyV2) + C^2yV2{QyVi)]' (5) 

The length, fleld, and time are, respectively, scaled by 
the coherence length ^o, ^c2,o = ^o/27r^o5 = 
47r^o ^^(t/c^ with the normal state conductivity crr^^'^^'^^ 
However, we here scale r]rn by r]o instead of r]o{T) = 
7^0 (1 — T/TcY^'^. r]Q is a uniform solution of r]i when 
V2 = 0, ^(r) = 1, and T = 0. The calculations are per- 
formed in a two-dimensional square area with each side of 
length of W and a current flow j, as schematically shown 
in Fig. n Outside the open boundary we set r]i =7/2 = 0, 
and at the lower and upper boundaries we set B{t) = 
and respectively, where H± = Ho^Wj/2 for an ap- 
plied fleld HqM 

We report the case when the gradient coupling is ab- 
sent {Ci2x = Ci2y = 0).2i We set Tc2/Tc = 0.9, = 4, 
C2 = 1, C3 = -C4 = 0.2, e = 0, and X = 0.5. If C4 is 
negative (i.e., a = tt), then ^1 = and 02 = ^^b, if the 
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FIG. 1: Configuration of the simulation. The current den- 
sity j flows in at the right-hand side boundary and out at 
the left-hand side boundary oi W x W square superconduc- 
tors. External fields and H- are applied outside the lower 
and upper boundaries, respectively. For the pinning case, we 
introduce pinning centers where the superconductivity is sup- 
pressed, represented by small solid squares. 



pairing is singlet. We note that our results for the vortex 
sheet do not significantly depend on the selection of the 
coeflficients in Eq. Q , because the appearance of the do- 
main wall only results from two states with the relative 
phase (a ± 7r)/2 being degenerate in free energy. We re- 
port results for T = O.lTc, H = 0.2i7c2,o, j = 2.6 x 10"^, 
and W = 57^0 as a typical example. In this simulation, 
we start from a uniform initial state with relative phase 
0, and switch on Hq and j at t = 0. 

In some cases of the penetration process of vortices, 
vortex sheet structures are created at the boundary and 
move toward the inside.^ On the other hand, the vor- 
tex sheet can be produced if we apply a magnetic fleld 
after preparing the domain wall at a zero fleld. In this 
case, some penetrating vortices are trapped at the do- 
main wall and these trapped vortices become a vortex 
sheet by splitting into half flux-quantum vortices. 

First, we study the free flux flow case without pin- 
ning centers, as shown in Fig. [21 The left panels show 
a color map of the spatial distribution of |7^i(r)| and 
1 7^2(1*) I . The black circle region shows a conventional vor- 
tex with a flux-quantum, where |7^i(r)| and 1 7^2(1*)! share 
the same vortex core. The green and red circle regions 
shows the r]i and 7^2 vortices, respectively, where only 
1 7^1 (r) I or 1 7^2(1*) I, respectively, has a vortex core and the 
other, 17^2(1*) I or |7^i(r)|, does not. These green and red 
vortex cores are located alternatively along a loop, form- 
ing a vortex sheet. The area including both an rji vortex 
and an 7^2 vortex is equal to that of a conventional vor- 
tex, indicating that each of the two vortices r]i and 7/2 
has half flux-quanta. The right panels shows the relative 
phase arg{7^i(r)/7^2(r)}, where we recognize that the vor- 
tex sheet is along the domain wall between the domains 
of relative phase tt (blue region) and relative phase (red 
region). The relative phase has a phase winding 27r or 
— 27r around the r]i or 7^2 vortex, respectively. In this free 
flux flow case, a conventional vortex flows upward, i.e., 
in a perpendicular direction to the current. The vortex 
sheet structure also simply flows with the same flow as 
conventional vortices. That is, the domain wall moves 
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FIG. 2: Time evolution of the vortex state in the free flux flow 
case without pinning centers. Snapshots of the total VK x VK 
area are presented for (a) t/54to = 17, (b) 22, and (c) 27. 
The left panels show color density plots of |?7i(r)| and |772(r)|. 
The black region is the core of the conventional vortex. The 
green and red regions are the cores of the 771 and 772 vortices, 
respectively, along the vortex sheet. The right panels show 
the relative phase arg{?7i(r)/?72(r)}. The red and blue regions 
are the domains of relative phase and tt, respectively. 



with the trapping half fiux-quantum vortices under cur- 
rent flow. The half flux-quantum vortices rarely escape 
from the domain wall of the vortex sheet in this free flux 
flow case. 

To study dynamical processes involving the vortex 
sheet and conventional vortices in a complicated man- 
ner, we introduce pinning centers to prevent the flux flow. 
Figure 01 shows the time evolution of the flux flow states 
in the presence of pinning, where pinning centers with 
a 1.4^0 X 1-4^0 square area are introduced periodically, 
as presented in Fig. Vortices move slowly, repeat- 
edly being trapped at pinning centers and escaping from 
them. The line of the vortex sheet also moves slowly and 
meanders due to the trapping of 771 and 772 vortices by 
the pinning centers. Compared with the slow motion of 
a conventional vortex and the line of the vortex sheet, 
r]i and r]2 vortices easily flow along the domain wall of 
the vortex sheet, as indicated by arrows in Fig. O^b). It 
seems that the barrier for escaping from a pinning center 




FIG. 3: Time evolution of the vortex state in the presence of 
pinning centers. Snapshots are presented for (a) t/54to = 196, 
(b) 276, and (c) 356. The left panels show color density plots 
of |?7i (r) I and |?72 (r) | . The right panels show the relative phase 
arg{7yi(r)/7y2(r)}. 



is lower along the line of the vortex sheet, because the 
vortex sheet is incommensurate with the pinning centers 
and the intervortex distance is short in this direction. 
Therefore, in the pinning case the domain wall forms a 
"channel" through which the vortex can flow. When the 
domain wall is connected to the sample boundary, vor- 
tices penetrate into the sample through the domain wall 
channel. These half flux-quantum vortices, quickly mov- 
ing through the channel, are emitted from the vortex 
sheet at the front corner of the slowly moving vortex sheet 
line. At the emitting point, pairs of half flux-quantum 
vortices combine to become a conventional vortex and the 
created conventional vortex is released from the domain 
wall of the vortex sheet. 

Figure S shows the time evolution of this emitting pro- 
cess in an enclosed region in Fig. O^b). In the upper panel 
in Fig. HJa) we see that vortices are trapped around pin- 
ning centers, represented by purple square regions. In the 
lower panels in Fig. ^ we find vortex centers of 771 or 772 
vortices as winding centers of the relative phase (between 
the yellow and purple regions). In Fig. EJa), r]i and 772 
vortices flow upward along the vortex sheet line. How- 
ever, the vortex sheet cannot easily move, because some 
of the r]i and r]2 vortices in the vortex sheet are trapped 
by pinning centers. Therefore, near the front corner of 
the domain wall shown in Fig. ^ the intervortex distance 
of r]i and r]2 vortices becomes short and the curvature of 
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FIG. 4: Time evolution of the creation of a conventional vortex from half flux-quantum vortices in the vortex sheet. Snapshots 
at (a) t/54to = 276, (b) 284, (c) 290, (d) 292, and (e) 296 are presented within the enclosed area shown in Fig. 3(b). The upper 
panels show color density plots of |?7i(r)| and |?72(r)|. The positions of the pinning centers are represented by purple squares in 
(a). The lower panels show the relative phase arg{?7i(r)/?72(r)}. The winding center between the yellow and purple regions is 
a vortex center of r/i or 772 vortices. 



the domain wall becomes sharp [(a)^(b)^(c)]. This de- 
formation of the domain wall shape helps the creation 
of the conventional vortex. After the core region of two 
772 vortices (the red regions in the upper panels) overlap, 
one of the 772 vortices forms a pair with a neighbor rji 
vortex [(c)^(d)]. This r]i and 772 vortex pair combines 
to become a conventional vortex (the black region) and 
leaves the domain wall [(e)]. We also observe the inverse 
process where a conventional vortex is trapped to the 
vortex sheet and changes to a pair of r]i and 7^2 vortices 
in the pinning case simulation. This event is rare, but 
occasionally observed. 

In summary, we have performed a TDGL simulation 
for the flux flow and the pinning in a two-component su- 



perconductor, and investigated the dynamics of the do- 
main wall and half flux-quantum vortices of the vortex 
sheet structure. The domain wall moves with the flux 
flow both in the free flux flow case and the pinning case. 
We succeeded in observing the creation process wherein 
a pair of half flux-quantum vortices is changed to a con- 
ventional vortex and released from the vortex sheet. Re- 
lating these phenomena, the vortex sheet structure may 
contribute to the flux flow resistivity or the magnetiza- 
tion process in multi-component superconductors, such 
as in PrOs4Sbi2, UPta, and Sr2Ru04. A quantitative es- 
timate of these interesting contributions is left for future 
studies. 
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